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We introduce an approach to studying a driven qubit-oscillator system in the ultrastrong coupling 
regime, where the ratio g/fl between coupling strength and oscillator frequency approaches unity or 
goes beyond, and simultaneously for driving strengths much bigger than the qubit energy splitting 
(extreme driving) . Both qubit-oscillator coupling and external driving lead to a dressing of the qubit 
tunneling matrix element of different nature: the former can be used to suppress selectively certain 
oscillator modes in the spectrum, while the latter can bring the qubit's dynamics to a standstill at 
short times (coherent destruction of tunneling) even in the case of ultrastrong coupling. 
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PACS numbers: 03.67.Lx, 42.50.Hz, 78.47.-p, 85.25.Cp 

The model of a two-level system coupled to a harmonic 
oscillator has been a standard applied in many different 
fields of physics. For instance, in quantum optics it is 
used to describe the interaction between light and mat- 
ter, leading to the field of cavity quantum electrodynam- 
ics (QED), where an atom interacts with the electromag- 
netic field of a resonator [l|, 0] . In the regime of strong 
coupling, where coherent exchange of excitations between 
atom and cavity is possible, those setups have become in- 
teresting for the field of quantum information with the 
atom being used as qubit and the cavity as information 
carrier. Additionally, the enormous progress in the field 
of circuit QED, where atom and cavity are replaced by 
superconducting circuits [3|45|, opens the door to the ul- 
trastrong coupling regime [6|-|8( with coupling strengths 
g between qubit and oscillator which are of the order of 
the oscillator frequency (typical values for cavity QED 
experiments are g/fl w 10 -6 ). Experimental realizations 
beyond the strong coupling regime have recently been 
reported 0, [l(J ■ The physics behind the qubit-oscillator 
system is usually analyzed within the Jaynes-Cummings 
model (JCM) which relies on a rotating- wave ap- 
proximation (RWA) with respect to g and provides deep 
insight into various effects of cavity and circuit QED. 
However, for the ultrastrong coupling regime the RWA 
fails, and theories beyond the JCM are needed, see, e.g., 
l_2j ■ An external probing of the qubit, e.g., by microwave 
radiation, can be modelled by the driven JCM 13[ , where 
a RWA is additionally invoked for the coupling between 
the qubit and the classical driving field, limiting the va- 
lidity of the model to moderate driving amplitudes. It 
has been shown 14[ that a strong external driving of the 
oscillator makes an inclusion of counter-rotating terms 
necessary even in the regime where the qubit's tunneling 
splitting equals the oscillator frequency (A = 0) and for 
couplings g/fl sw 0.1, parameters for which the JCM is 
commonly used in the undriven case. Similar effects are 
expected if instead of the cavity the atom is driven. Such 
extreme driving strengths have already been experimen- 
tally realized (l5l - fl8| leading to a dressed qubit state [lj| . 
In this work we examine analytically the spectrum and 



dynamics of a system exposed to both ultrastrong cou- 
pling and extreme driving. To go beyond common RWA 
schemes, the qubit is treated within Floquet theory jiol ). 
while the coupling to the quantized field is included using 
a polaron transformation [2l|. We study various control 
possibilities of the qubit's dynamics, from tunneling sup- 
pression to selective frequency generation. 
The spectrum. The Hamiltonian of the driven qubit- 
oscillator system reads 

H = -(h/2)[e(t)a z +Aa x ]+hg& z {B^+B)+hriB J 'B, (1) 

where &i are the Pauli matrices, and B\ B the cre- 
ation and annihilation operators of the oscillator. A si- 
nusoidal variation of the static bias e is given by e(t) = 
e + A cos tu cx t. In order to treat the time dependence, we 
analyze the Floquet Hamiltonian T-L — H(t) — ihdt and 
consider the extended Hilbert or Sambe space fj (g> T [2(J , 
where Sj is the Hilbert space of the undriven system and 
T the space of the time-periodic functions. A basis of 1 is 
provided by the vectors |/), whereby (t\l) = expj— ilui K yt\. 
The eigenstates of the driven qubit for A = are [22| : 
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with quasienergies fe9 



function J n (x) 
the oscillator, 
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. For no interaction between the qubit and 
g = 0, this eigenbasis is easily extended 



to the full Hamiltonian ri by n K }} = n ))\K), 
with \K) being an eigenstate of the oscillator. For the 
coupled system (g ^ 0), the eigenstates can be found 
for A = with the help of the polaron transformation 
U = cxp{.g(S - B^)a z /n} and are a combination of the 
Floquet states of the qubit and displaced oscillator states: 
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|t /i,n t K)) = cxp [±[g(B - 
with the quasienergies 

- T he/2 - ftnu cx + hKfl - hg 2 /fl. (4) 
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For A = 0, this result is analytically exact and treats the 
problem for arbitrary coupling strength g. Figure [T] shows 



At e 
numbers K 



the energy spectrum of Eq. (|4]) for A = (blue squares) . 

mWgx — Lfl crossings occur; the quasienergies 
and K are degenerate, with the integer 
+ L = 0, 1, . . . , oo and m,n = — oo, . . . , oo 
[2^ | . Note that for L ^ there are always L nondegen- 
erate levels. For L > those are the first L spin-down 
states (positive slope), while for L < the first L spin- up 
states (negative slope) . At finite A avoided crossings oc- 
cur in the energy spectrum at the sites of the resonances 
(red triangles and black dots in Fig. [lj. 
To explain the origin of these avoided crossings we ex- 
press H in the basis (O yielding the off-diagonal elements 
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,K 



{{l,n,K\Aa x \t,n',K>)) 
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The dressing A m = AJ m (A/u) ex ) of the tunneling ma- 
trix element results from the external driving [20I ] , while 
a L/2 y/K\/(K + L)\ L^\a)e-^ stems from 



the coupling to the oscillator [23M25I ] with L^'{x) be- 
ing the K th generalized Lagucrre polynomial and a = 
(2.9/tt) 2 . 

In order to calculate the energy spectrum for finite A, 
we make use of Van Vleck perturbation theory in anal- 
ogy to 22 1 and [25j]. To first order in A, we only take 
into account states degenerate for A = in 7i, together 
with the matrix elements connecting them. Corrections 
from the remaining off-diagonal elements are calculated 
to second order. The resulting effective Hamiltonian con- 
sists of 2 x 2 blocks (without loss of generality we assume 
that L > 0): 
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where we introduced the second-order corrections 



/(e+puj cx ±Pn). 



(2) = v V (K n ' K 
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{p,F}/{-m,±L} 

(7) 

The dressed tunneling matrix element in Eq. © deter- 
mines to first order the width of the avoided crossings 
in Fig. [TJ Dominant crossings are found for e = mw ex , 
where the static bias is an integer multiple of the driving 
frequency, and thus L = 0. That means that both states 
belong to the same oscillator quantum number K, and 
the dressing contains a Laguerre polynomial of the kind 
L° K (a). For the 2x2 block in Eq. © the eigenvalues to 
the eigenstates \^^ n £ K )) are found easily: 
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FIG. 1: (Color online) Quasienergy spectrum of the qubit- 
oscillator system against the static bias e for weak coupling 
g/ojox = 0.05. Further parameters are A/tj ex = 0.2, Q./ui cyi = 
\/2, A/uj cx = 2.0. The first six oscillator states are included. 
Numerical calculations are shown by red (light gray) triangles, 
analytical results in the region of avoided crossings by black 
dots. A good agreement between analytics and numerics is 
found. Blue (dark gray) squares represent the case A = 0. 



where the upper indices denote the state of the qubit, 
the Floquet mode and oscillator quantum number, while 
the lower indices stand for the resonance condition. The 
width of the avoided crossings is given by 



Q n ' K r = 

m.L 



e - mu cx + m+ (ef^ K+L + s\% >K )/4 
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the dressed oscillation frequency, which, together with 
Eq. (|5J , is one major result of this work. For the L non- 
degenerate spin-down states the quasienergies and eigen- 
states are simply E^ K + \s^ n K and ||, n, K)). 
FigureOshows the quasienergy spectrum against the cou- 
pling strength g. For simplicity, we study the unbiased 
case e = 0, which implies m = L = and hence gaps 
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Thus, for g 







and A ^ 0, the twofold degeneracy of the unperturbed 
case is lifted by a gap of width Ao. For j ^ 0, the gap 
size is further determined by the Laguerre polynomial, 
so that additional degeneracies can occur at the zeros of 
L K (a). When choosing the driving amplitude A such 
that Ao = the twofold degeneracy is kept for arbitrary 
g and K. Because the dressing by the Bessel function 
does not depend on g or the oscillator level, we reach the 
remarkable conclusion that the coherent destruction of 
tunneling (CDT), predicted for a driven qubit |26j |. might 
occur also for a qubit-oscillator system in the ultrastrong 
coupling limit. In Fig. 02 the dressed oscillation frequen- 
cies are plotted against the dimensionless coupling g/fl. 
Next to an exponential decay, they exhibit zeros that de- 
pend through the Laguerre polynomial characteristically 
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FIG. 2: (Color online) Quasienergy spectrum against the cou- 
pling strength g in the unbiased case e — 0. Further, we set 
u) ex /Q = 5.3 and A/Q = 1.0. The spectrum is examined 
for A/Q = 8.0 and A/Q = 12.74 . For the former, avoided 
crossings of amplitude K occur, which vanish at values of g 
yielding zeros of the Laguerre polynomials. For the latter, all 
Q vanish simultaneously for all values of g, since the CDT 
condition Jo(A/(jj cx ) — is fullfilled independently of the cou- 
pling strength g. As a reference the case A = is shown. 

on the oscillator quantum number K. Hence, because 
the qubit's dynamics involves several oscillator levels, we 
predict that suppression of tunneling cannot be reached 
by just tuning the coupling g. 

The dynamics. To prove the statements above, we cal- 
culate the survival probability of the qubit Pj,_>j,(t) := 
(X \{>red{t)\ 4), where /3 rc d is obtained by tracing out the 
oscillator degrees of freedom from the density operator 
of the qubit-oscillator system: 

/Cr K V>,o do) 

wil li j, ■■= E a m K L ~E? m % and {a, /?} e {-, +4}- 
The index n has been dropped, because it just leads to 
an overall phase and thus has no influence on the dynam- 
ics. The time-dependent Floquet modes |f^(t)) can be 

determined from the eigenstates | < J > ^ 1 ' i £)) of the effective 
Hamiltonian ^ (22[. In Fig. |H we plot the dynamics for 
zero bias and three different coupling strengths, indicated 
by the dashed lines in Fig. [3] We assume factorized start- 
ing conditions for t = 0, with the qubit in the state | J.}, 
and the oscillator in thermal equilibrium obeying a Boltz- 
mann distribution. From Eq. (|10j) , one expects two main 
oscillatory contributions, namely, u>^ K £ K = ±0^ L and 

LO 'rnL lK = {K — K')Q. Also sums of both can occur. For 
weak coupling g/Q = 0.1 (a), the analytical calculation 
shows oscillations between the states | \) and | f) with 
the single frequency Q°. For stronger coupling g/Q = 0.5 
(b), a second small peak at Q 2 occurs in the Fourier spec- 
trum, whose effect on the survival probability is almost 




FIG. 3: (Color online) Size of the avoided crossing K against 
the dimensionless coupling strength g/Q for an unbiased qubit 
(e = 0). Further, A/Q = 0.4, cj cx /fi = 5.3 and A/Q = 8.0. 
Q K vanishes at the zeros of the Laguerre polynomial L° K (a). 
The dashed lines (a), (b), (c) represent g/Q = 0.1, 0.5, 1.0, 
respectively, as considered in Fig. [3] 
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FIG. 4: (Color online) Dynamics of the qubit for e — 0, 
A/51 = 0.4, u eK /Q = 5.3, A/Q — 8.0, and temperature 
hQ(ksT)^ 1 = 10. The graphs show the Fourier transform 
F(v) of the survival probability Pj._^^(t) (see the insets). We 
study the different coupling strengths indicated in Fig. [3] 
g/Q — 0.1 (a), 0.5 (b) and 1.0 (c). Analytical results are 
shown by black curves, numerics by dashed orange curves. 

not visible. The peak at fi 1 is absent, because the cor- 
responding Laguerre polynomial vanishes at this value 
exactly, see Fig. [31 Correspondingly, we observe a tun- 
neling reduction compared to case (a). In Fig. |4]Jc) we 
are with g/Q, = 1.0 already deep in the ultrastrong cou- 
pling regime. The frequency Q 1 is now different from 
zero, and additionally f2 3 appears. The lowest peak be- 
longs to the frequencies O 2 , and f2 4 , which are equal 
for g/Q = 1.0, see Fig. [3j A complete population in- 
version again takes place. Our results are confirmed by 
numerical calculations. For g — 0.5, 1.0, the latter yield 
additionally fast oscillations with Q and w GX - Further- 
more, Q 1 is shifted in Fig. @|c) slightly to the left, so 
that concerning the survival probability the analytical 
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FIG. 5: (Color online) Coherent destruction of tunneling in 
a driven qubit-oscillator system. The same parameters as 
in Fig. [4] are used except that A/Q — 12.7, which leads to 
Ao = 0. Three coupling strengths are examined: g/Q = 0.1 
(a), 0.5 (b) and 1.0 (c). The analytical calculations (black, 
dashed lines) predict complete localization for all three cases. 
Also the numerics (red curves) shows strong localization for 
short timescales with fast oscillations overlaid. For long times 
this localization vanishes (see inset in (a)). 



and numerical curves get out of phase for longer times. 
To include also the oscillations induced by the driving 
and the coupling to the quantized modes, connections 
between the degenerate subspaces need to be included in 
the calculation of the eigenstates of the full Hamiltonian 



22, 



While tuning the coupling g to a zero of a Laguerre 
polynomial corresponding to a dominant oscillator mode 
yields a reduction of tunneling, tuning the driving am- 
plitude A to a zero of a Bessel function can yield almost 
complete localization at short times. As already noticed 
in Fig. [21 this phenomenon is independent of the coupling 
strength g. We choose in Fig. [5]the driving amplitude A, 
so that Ao = 0. This is the same condition as found for 
CDT in a driven qubit [261 ], Analogously, our analytical 
solution now predicts localization for arbitrary coupling 
strength g. All dressed oscillation frequencies Vl K van- 
ish. However, third-order corrections in A will give small 
contributions to Ao 27(. Hence, a numerical exact so- 
lution yields oscillations of Pj-^t) with a long period. 
On a short timescalc and for oj cx ^ A also the numerical 
solution appears to be strongly localized, while for long 
times, the inset in Fig. [5] (a) shows complete population 
inversion for the numerics. 

In conclusions, we developed a powerful formalism to 
investigate analytically a qubit-oscillator system in the 
ultrastrong coupling and extreme driving regime, a sit- 
uation which is in close experimental reach and offers 
excellent control possibilities. Our approach relies on 
perturbation theory with respect to a single parameter 
only, the qubit tunneling matrix element A, and thus 



goes beyond the driven Jaynes-Cummings model, with 
no rotating-wave approximation being applied. 
We acknowledge financial support under DFG Program 
SFB631. We thank Sigmund Kohler for helpful remarks. 
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